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The dynamical heterogeneity in supercooled liquids measured by a molecular dynamics simulation 
has been quantified on the basis of the multifractal formalism. The singularity spectrum becomes 
broader as the glass transition is approached. This behavior is similar to that observed in the case 
of the Anderson-localization transition. 
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The dynamical heterogeneity in supercooled liquids has been studied intensively in recent years. It can be quantified 
by introducing the concept of bond-breakage and the distribution of the broken bonds near the glass-transition is found 
to be similar to that of the spin state of the Ising model near the critical point. The bond-breakage corresponds 
to the collective jump motion of particles among cages. 

On the other hand, we have analyzed the dynamical heterogeneity of the critical spin state of the Ising model 
on the basis of the multifractal formalism. || In this Short Note we try to combine above two approaches and 
quantify the dynamical heterogeneity observed in a molecular dynamics simulation of a supercooled liquid in terms of 
multifractality. The so-called cooperatively rearranging regions (CRR) can be identified with weakly bonded regions 
with various sizes which form multifractal structure. 

We use the same model as in Ref. 1, the binary mixture of soft core particles. Two species of particles, 1 and 2, 
with the size <j\ and 02, and the mass mi and mi are considered. The interaction potential is given by 

v aP {r)=e-C-ff\ (1) 

where r is the distance between two particles and a a /3 — (a a + <Jp)/2 with a, f3 = 1, 2. 

The molecular dynamics simulation is done using the second order symplectic integrator [Q and the microcanonical 
ensemble. The temperatures in the following are determined by the kinetic energy. 
The bond is introduced as in Ref. 1 when 

rij(to) < 1-1^/3, (2) 

where rij(to) is the distance between i-th particle of species a and j-th particle of species (3 at time to- The position 
of the bond is represented by the midpoint of the two particles at this time. The bonds are unambiguously defined 
0j owing to the sharpness of the first peak of their pair correlation function. The bond is regarded as surviving at 
time t + 1 when 

rij(to + t) < 1.6a af3 . (3) 

The number of surviving bonds N suv (t + t) is fitted by 

N suv (t Q +t) = N suv (t ) exp(-i/r b ), (4) 

introducing the life-time of bonds t^. In the following the broken bonds are determined for t = 0.05rb. 

In our simulation the followings are fixed. Two species of particles have the same particle number, Ni = N2 = 80000. 
The size and mass ratios are chosen as (J%la\ = 1.4 and mijm\ — 2 respectively. The time step of the integration 
At is chosen as At = 0.005ti with t\ = (miaf/e) 1 ^ 2 . The potential is truncated at r — 4.5<ti. The particle density is 
chosen as n — 0.8 /a 2 . The simulations were performed on our personal computer with Pentfum-4 CPU. 

In Fig. 1 the snapshot of the broken bonds is shown where the position of the bond is defined at time to as mentioned 
in the above. In Figs. 1(a) and 1(b) the temperatures are chosen as T = 2.54 and T = 0.859 respectively in the unit of 
e/fce- These two states are realized starting from an initial state with the temperature 9.69 and cooled in a stepwise 
manner as in Ref. 1. In the stepwise cooling process the temperatures in the first, second, third and fourth stages 
are 9.69, 2.54, 1.43 and 0.859 respectively. For the former state the temperature was changed from 9.69 to 2.54 after 
spending 5000Ai at the first stage and the broken bonds are observed at the second stage with to = 5000Ai. For the 
latter state the temperature was changed successively from 9.69, 2.54, 1.43 to 0.859 after spending 5000Ai at the first, 
second and third stages each and the broken bonds are observed at the fourth stage with t = 60000At. The times 
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to and t are measured from the instance of the temperature change to the stage of the observation. The life-times of 
bonds were measured to be r b = 4.09 x 10 3 Ai for T = 2.54 and r b = 2.43 x 10 5 Ai for T = 0.859. 

Qualitatively the snapshot shows a multifractal behavior and the degree of intermittency increases as the temper- 
ature is decreased. In the following we quantify this observation on the basis of the multifractal formalism. 

The multifractality is quantified by the singularity spectrum, /(a), and the spectrum is obtained by the following 
formulae, || 

a(q) =J2 Mfc(9,0k/i fc (l, 1)1 HI/ L), (5) 

k=l 

and 

Nl 

f(a(q)) = J2^(q,l)^^ k (q,l)/hx(l/L), (6) 

k=l 

with 

Nl 

w(?.!) = {w(i)}7E(w( ! )}' ( ? ) 

k' = l 

Here the measure, fJ.k(l), is defined as the probability to find a broken bond in the fc-th square box of the linear 
dimension I in a snapshot of broken bonds. The total number of the boxes Nl = (L/l) x (L/l) where L is the linear 
dimension of the simulation region with periodic boundary condition. 

In Fig. 2 the singurarity spectra for T = 2.54 and T — 0.859 are shown. The dots represent the data set 
(oi(q), f(a(q))) determined by the above formulae and should be tangential to f(a) = 2 and f(a) — a. The maximum 
a max and minimum a min are determined by the intercept of the spectrum extrapolated to the a-axis and correspond 
to the singularity exponent for the box with the lowest density of broken bonds and the one with the highest density 
respectively. 

From the data in Fig. 2 we can quantify the degree of intermittency. The singularity spectrum becomes broader 
as the temperature is decreased approaching the glass transition. At the same time a max becomes larger and a m ; n 
becomes smaller. Namely a max = 2.74 and a m j n = 1-66 for T = 2.54 and a max — 2.97 and a m i„ = 1.57 for T = 0.859. 
This finding is similar to that reported in the case of the Anderson localization where the singularity spectrum 
becomes broader as the strength of the impurity potential is increased and the localization transition is approached. 
In both cases the degree of intermittency corresponds to that of localization. 

In conclusion we have quantified the dynamical heterogeneity in a model supercooled liquid measured by a molecular 
dynamics simulation on the basis of the multifractal formalism. In supercooled liquids various sizes of CRR are formed 
and constitute multifractal structures. The degree of intermittency increases approaching the glass transition. Since 
our simulation parameters are still far from those for the glass transition, we need further study to elucidate the 
nature of the glass transition. 
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Science and Technology of Japan. 



[1] This paper has been submitted to the Short Note section of J. Phys. Soc. Jpn. 

[2] R. Yamamoto and A. Onuki: J. Phys. Soc. Jpn. 66 (1997) 2545. 

[3] W. Sakikawa and O. Narikiyo: J. Phys. Soc. Jpn. 71 (2002) 1200. 

[4] M. Tuckerman, B. J. Berne and G. J. Martyna: J. Chem. Phys. 97 (1992) 1990. 

[5] M. Schreiber: Computational Physics, ed. K. H Hoffmann and M. Schreiber (Springer- Verlag, Berlin, 1996) p. 147. 



FIG. 1. The snapshot of the broken bonds for (a) T = 2.54 and (b) T = 0.859. 
FIG. 2. The singularity spectrum for T = 2.54 and T = 0.859 with L/l = 30. 
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